Spectroscopy on a single trapped 137 Ba + ion 
for nuclear magnetic octupole moment determination 
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We present precision measurements of the hyperfine intervals in the 5D 3 / 2 manifold of a single 
trapped Barium ion, 137 Ba + . Measurements of the hyperfine intervals are made between tuf = 
sublevels over a range of magnetic fields allowing us to interpolate to the zero field values with an 
accuracy below a few Hz, an improvement on previous measurements by three orders of magnitude. 
Our results, in conjunction with theoretical calculations, provide a 30-fold reduction in the uncer- 
tainty of the magnetic dipole (A) and electric quadrupole (B) hyperfine constants. In addition, we 
obtain the magnetic octupole constant (C) with an accuracy below 0.1 Hz. This gives a subsequent 
determination of the nuclear magnetic octupole moment, Q, with an uncertainty of 1% limited al- 
most completely by the accuracy of theoretical calculations. This constitutes the first observation 
of the octupole moment in 137 Ba + and the most accurately determined octupole moment to date. 

PACS numbers: 32.10.Fn, 21.10.Ky, 27.60.+j 



High precision measurements of the hyperfine structure 
provides stringent tests for state-of-the-art atomic struc- 
ture calculations. These calculations play a crucial role in 
the interpretation of parity nonconservation (PNC) ex- 
periments which provide important tests of the standard 
model at low energy [1] . The accuracy of calculated PNC 
matrix elements can be assessed by comparing measured 
hyperfine structure constants with calculated values [2J. 
In addition, the hyperfine structure provides insight into 
the nuclear structure of atoms [3]. 

In this paper we present precision measurements of 
the hyperfine intervals of the 5D 3 / 2 manifold of 137 Ba + . 
By combining high precision radio frequency (rf) spec- 
troscopy with shelving techniques [11 [5] on singly trapped 
ions, we measure the hyperfine intervals of the 5D 3 / 2 
manifold to an accuracy below a few Hz. Together with 
theoretical calculations, this permits a 30-fold reduction 
in the uncertainty of the magnetic dipole (A) and elec- 
tric quadrupole (B) hyperfine constants and the first ob- 
servation of the magnetic octupole moment in 137 Ba + . 
We determine the magnetic octupole moment, Q, to an 
accuracy of 1% percent limited almost entirely by uncer- 
tainties in the theory. 

Ba + is an excellent candidate for magnetic octupole 
determination. Long lived metastable D states permit 
high precision spectroscopy measurements of the hyper- 
fine levels, and high precision calculations are possi- 
ble [BJ. In [7] it was shown that the hyperfine intervals, 
SWf — Wf — Wf+i, of the 5D 3 / 2 manifold can be writ- 
ten 

SWo = -A + B-56C-^r ) + ^-, (1) 
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SWi = -2A + B + 28C - —77, (2) 



SW^-SA-B-SC-^r,-^ (3) 

where C is the magnetic octupole hyperfine constant 
while r\ and £ are the correction terms characterizing the 
mixing with the upper 5D 5 / 2 manifold. With measured 
values of the hyperfine intervals, SWk, and theoretical 
estimates for the correction terms we can solve the equa- 
tions for the hyperfine constants. The octupole moment, 
f2, can then be extracted from our theoretical result [8] 

C = -0.581(6) [ — ^— ) kHz, (4) 
\fi N x b/ 

using the CCSD(T) method described in [2J[5]> where /in 
is the Bohr magneton and b is the barn unit of area. 

The procedure to measure the hyperfine intervals is 
similar to that proposed in [5J UOj and the relevant level 
structure is given in Fig. [T] The ion is first Doppler 
cooled and then optically pumped into the state \F = 
2, nip — 2) of the 6Si/ 2 level. To measure the hyperfine 
interval, SWk, the ion is then shelved to the state \F" = 
k,mp =0) of the 5D 3 / 2 level using a two photon Raman 
transition similar to the one used in [5]. An rf antenna 
is then turned on to drive F" -R> F" + 1 transitions. The 
signal generator used for the rf antenna is synchronized 
with a GPS-disciplined Rb clock. 

To determine if the hyperfine transition occurred we 
use a second Raman pulse to transfer the state \F" — 
k,m,F = 0) back to the \F — 2,mp — 2) ground state 
and then optically pump on the 6S!/ 2 to 6P 3 / 2 transi- 
tion at 455 nm. If the hyperfine transition occurred, the 
ion will remain in the \F" = k + 1, mjr = 0) state of the 
5D 3 / 2 level, otherwise it will be shelved to the 5D 5 / 2 level 
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6Si/ 2 

FIG. 1. Relevant levels of 137 Ba + for the rf spectroscopy: 
The ion is prepared in the \F = 2, rriF = 2) state of the 6S1/2 
level from where it is shelved to \F" ,mpn = 0} of the 5D 3 / 2 
level with a pair of Raman beams which are red-detuned by 
« A = 2tt x 500 GHz from the 6P1/2 level. The rf transition 
is detected by shelving to the 5D 5 / 2 level (see text). 



with sa 88% efficiency. Subsequent driving of the 6S1/2 to 
6P1/2 an d 5D3/2 to 6P!/ 2 transitions using the Doppler 
cooling beams then provides a fluorescence measurement 
for the probability of driving the hyperfine transition: the 
ion being bright if the hyperfine transition took place and 
dark otherwise. The detection efficiency of this scheme 
is limited to approximately 88% by the branching ratio 
between the 6P3/2 and 5D 3 / 2 which results in unwanted 
population of the 5D 3 / 2 level when optically pumping to 
the 5D 5 /2 level. This does not impact on the accuracy 
at which we can measure the hyperfine transition prob- 
ability, but only on the amount of averaging needed to 
achieve a particular level of accuracy. 

The experiments are performed in a four-rod linear 
Paul trap, similar to the ones described in [SI [TT]. The 
trap consists of four stainless steel rods of diameter 
0.45 mm whose centers are arranged on the vertices of 
a square with 2 mm length of the side. A 5.3 MHz 
rf potential with an amplitude of 170 V is applied via 
a step-up transformer to two diagonally opposing elec- 
trodes. A small DC voltage applied to the other two 
electrodes ensures a splitting of the transverse trapping 
frequencies and rotates the principle axes of the trap with 
respect to the propagation direction of the cooling lasers. 
Axial confinement is provided by two axial needle elec- 
trodes separated by 2.4 mm and held at 33 V. Using 
this configuration, the measured trapping frequencies are 
{lo x , u y ,uj z )/2Tr « (1.7, 1.5, 0.5) MHz. 

Doppler cooling and detection is achieved by driving 
the 6S1/2 to 6P!/ 2 transitions at 493 nm and repump- 
ing on the 5D 3 / 2 to 6P1/2 transitions at 650 nm. The 
493 nm laser is passed through two electro-optic modu- 
lators in order to generate the sidebands required to ad- 
dress all possible transitions between the 6S1/2 to 6P1/2 



levels. Similarly, to address all the 5D 3 / 2 to 6P1/2 tran- 
sitions, the repumping laser at 650 nm is split into four, 
frequency shifted by acousto-optic modulators, and then 
recombined into a single fiber. All laser fields are linearly 
polarized perpendicular to a magnetic field of approx- 
imately 1 G. This configuration avoids unwanted dark 
states in both the cooling and detection cycles. Addi- 
tionally, the field ensures a well defined quantization axis 
for optical pumping and state preparation. 

The magnetic field also gives a second order Zeeman 
shift of the \F",mp = 0) levels which must be ac- 
counted for when measuring the hyperfine intervals. We 
achieve this by additionally measuring transitions be- 
tween \F",m F = 0) and \F" + l,m F = ±1). Half of 
the difference frequency between the Amp = +1 and 
the ArriF = —1 transitions is then given by the lin- 
ear shift ^B9FB/h with any quadratic shifts cancelled. 
By measuring the Amp = transition over a range of 
magnetic fields, each calibrated by measurements of the 
Amp — ±1 transitions, we can map out the full field de- 
pendence of the Amp — transition and interpolate to 
the desired zero field result. This also yields the second 
order Zeeman shift coefficient for the Amp = tran- 
sitions providing a useful consistency check within our 
measurements. 
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FIG. 2. Transition probability as a function of the external 
rf frequency for the \F" = 0,m F „ = 0}-<^\F" = l,m F » = 0) 
transition. 

An example of a hyperfine splitting measurement for 
a fixed value of the magnetic field is shown in Fig. [2] 
where the transition probability for \F" = 0, mp = 0) to 
\F" — 1, mp = 0) is plotted as a function of the external 
rf frequency. The rf drive power is adjusted to give an 
on resonant 7r-pulse time of approximately 10 ms result- 
ing in a resonance width (full width half maximum) of 
approximately 50 Hz. The resonance is scanned in 2 Hz 
frequency steps, and the probability of undergoing the 
hyperfine transition is determined by the average over 
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200 measurements. The data is fitted [H] via a \ 2 min- 
imization to a sine function with additional offset and 
amplitude parameters to account for imperfect shelving. 
From the fits we also extract the 68% confidence limits 
on the fit parameters giving a determination of the center 
transition frequency with an accuracy of about 1 Hz. At 
this same fixed magnetic field, a similar measurement is 
also performed for the Amp = ±1 transitions. For this 
case we use a 7r-pulse time of 0.5 ms and scan the reso- 
nances in 20 Hz steps. This gives a measurement of the 
resonant frequency for these transitions with an accuracy 
of about 20 Hz which corresponds to a field accuracy of 
about 50 /iG. 

The time to perform the Amp = and field calibra- 
tion measurements at a particular value of the magnetic 
field is about 15 min. It is therefore necessary to con- 
sider the effects of magnetic field drifts which give rise to 
additional errors in the field calibration. Hence we have 
monitored the variation of the magnetic field over the 
course of a day by measuring the first order Zeeman shift 
of the Amp = +1 transition. From these measurements 
we estimate the one standard deviation error in the field 
calibration to be approximately 450 fiG which dominates 
the 50 /iG error extracted from the Amp — ±1 field cal- 
ibration measurements. 

Measurements for the hypcrfine interval, SWq, are 
shown in Fig. [3] For each hyperfine interval we have 
taken two sets of data on separate days. For each set, 
we measured the transition frequency at 10 values of the 
field: 5 with the field above zero and 5 below. The in- 
set highlights two points taken at a similar field setting 
on separate days. Confidence limits from the fits to the 
resonance scans determine the vertical errors, which are 
smaller than the thickness of the lines shown in the inset. 
We fit the data for each hyperfine interval to a quadratic 
form, a^B 2 — SWk, using a x 2 minimization. Since the 
second order Zeeman shift coefficients, ak are determined 
by the three values of SWk, we fit all three quadratic 
forms simultaneously. However, we may still determine 
a x 2 statistic for each data set. Since the ak are only 
weakly dependent on the SWk, the minimization proce- 
dure is equivalent to three independent single parameter 
fits with the ak fixed to values consistent with the fitted 
values of SWk- In Table [I] we give the SWk along with the 
reduced \ 2 for each fit. Errors reported here are again 
the 68% confidence limits extracted from the fits. For the 
field calibration and calculation of the ak we have used 
values of gj and gj reported in [13] and [14] respectively. 

It is worth noting that the quadratic form used for the 
fitting is only an approximation. As the field strength 
increases, higher order terms become important, which 
can shift the zero field value extracted from a quadratic 
fit. The smaller the hyperfine splitting, the larger the 
effect. It is for this reason we have restricted our field 
values to ~ 1 G. Over the range of field values we have 
used, we estimate that this effect could be as large as ~ 
0.5 Hz for SWq. For field values in the range 1 ~ 2 G this 
systematic error would be comparable to the uncertainty 
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FIG. 3. Measured hyperfine splitting as a function of the 
magnetic field for the \F" = 0, m F » = 0}o|F" = 1, m F » = 0) 
transition. 



in our current measurement . 

Additional systematic errors could arise from three 
possible sources: off-resonant coupling of the rf drive to 
other Zeeman levels, coupling of the ion trap rf field to 
the ion via micromotion, or AC Stark shifts due to leak- 
age light from the 650 nm lasers used in the experiment. 
In Table [H] we give estimates for these shifts. For the 
rf drive field, the largest effect comes from the field cali- 
bration measurements since the Rabi rate for these mea- 
surements is significantly larger than for the Amp = 
measurements. In this case, the Stark shifts vary with 
the magnetic field and can alter the quadratic form in 
Fig. [3] causing a shift in the zero- field point. This effect 
is largest for the F = O F = 1 transition and the esti- 
mate given in the table assumes equal power available for 
the Amp = 0, ±1 transitions. The effects of micromotion 
were estimated using the results in [TS] . Additionally, we 
shifted the ion position to induce micromotion and saw 
no measurable shift in the resonance. For the 650 nm 
shelving beam, it was found that 1 fiW of leakage light 
was enough to shift the resonances by ~ 60 Hz consis- 
tent with calculations. The estimate given in Table pi] is 
based on this measurement and the additional attenua- 
tion used to remove the shift. For the repumping beams 
we attempted to measure the leakage light level by mea- 
suring the rate at which population was removed from 
the D 3 / 2 level. As there was no statistically significant 
population loss over 200 ms, the estimate given is an up- 
per bound. 

The hyperfine coupling constants inferred from the 
measured values of SWk are given in table |III| To account 
for the possible systematic errors we have added 0.5 Hz 
to the errors given in table [I] The correction terms were 
calculated from the expressions given in [S] , using the ma- 
trix elements (5Z? 3/2 1 |Ti | |5Z? 5/2 > = -995(10) MHz//i N 
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and (5Z? 3/2 ||r 2 ||5L> 5/ 2) = 339(5) MHz/b 0. The mag- 
netic moment 0.937365(20)/in and quadrupole moments 
246(l)b where taken from [Tf)] and [S] respectively. 

The final values obtained for A and B are within three 
standard deviations of the results reported in |17j with 
a 30-fold reduction in the uncertainty. In addition, our 
value for C has an accuracy of 0.2% and is the first ob- 
servation of the magnetic octupole moment in 137 Ba + . 
Using Eq. Q we calculate an octupole moment of 

r2( 137 Ba+) = -0.06290(67) (/a n x b). 

This, to our knowledge, is the most accurately deter- 
mined octupole moment to date. 



TABLE I. Measured hyperfine intervals, SWk, for the 5D 3 / 2 
manifold of 137 Ba+. 



TABLE III. Hyperfine coupling constants. 



Transition 


—5W (Hz) Reduced X 2 


F=0 -» F=l 


145 193 549.3 (2.8) 1.01 


F=l F=2 


334 921 347.13 (89) 1.60 


F=2 -» F=3 


613 730 628.08 (22) 0.47 


TABLE II. 


Estimates of the possible systematic errors. 


Source 


Error estimation 



Off resonant rf coupling 
Micromotion 
Stray Raman light 
Stray repump light 
Higher order terms (SWo) 



~ 0.1 Hz 
~ 5 mHz 
~ 60 mHz 
< 0.5 Hz 
~ 0.5Hz 



A (Hz) 



B (Hz 



C (Hz) 



Uncorr. 189730524.90(31) 
77 corr. 805(16) 
C corr. -228.5(4.2) 
Corr. 189731101(17) 



44538793.6(1.0) 
-1610(32) 
-571(10) 
44536612(34) 



32.465(42) 

4.081(75) 
36.546(86) 



In conclusion we have performed high precision mea- 
surements of the hyperfine structure in 137 Ba + . Our 
measurements have greatly reduced the uncertainty of 
the currently available hyperfine structure constants and 
have provided an estimate of the nuclear octupole mo- 
ment, f2, accurate to 1%. We note that our inferred 
value of has the opposite sign compared to the esti- 
mate given in 7 . This estimate is based on the shell 
model for which such discrepancies have also been re- 
ported for Cesium [18] and Rubidium [19]. Thus our 
measurements may present an interesting test case for 
comparison with nuclear-structure calculations. Approx- 
imately 10% of the estimated value for 17 is given by the 
correction factor, £. Future planned measurements on 
the 5D 5 /2 manifold of the same ion will allow the depen- 
dence on this factor to be eliminated. This would also 
provide additional consistency checks between measure- 
ments and calculations. 
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